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ABSTRACT 

Using metallicities from the literature, combined with the Revised Bologna Catalogue of pho- 
tometric data for M3 1 clusters and cluster candidates (the latter of which is the most com- 
prehensive catalogue of M31 clusters currently available, including 337 confirmed globular 
clusters - GCs - and 688 GC candidates), we determine 443 reddening values and intrinsic 
colours, and 209 metallicities for individual clusters without spectroscopic observations. This, 
the largest sample of M3 1 GCs presently available, is then used to analyse the metallicity dis- 
tribution of M31 GCs, which is bimodal with peaks at [Fe/H] « —1.7 and —0.7 dex. An 
exploration of metallicities as a function of radius from the M3 1 centre shows a metallicity 
gradient for the metal-poor GCs, but no such gradient for the metal-rich GCs. Our results show 
that the metal-rich clusters appear as a centrally concentrated spatial distribution; however, the 
metal-poor clusters tend to be less spatially concentrated. There is no correlation between lu- 
minosity and metallicity among the M3 1 sample clusters, which indicates that self-enrichment 
is indeed unimportant for cluster formation in M3 1 . 

The reddening distribution shows that slightly more than half of the GCs are affected by a 
reddening of E(B - V) < 0.2 mag; the mean reddening value is E(B - V) = 0.28^4 
mag. The spatial distribution of the reddening values indicates that the reddening on the north- 
western side of the M3 1 disc is more significant than that on the southeastern side, which is 
consistent with the conclusion that the northwestern side in nearer to us. 

Key words: galaxies: individual (M31) - galaxies: star clusters - globular clusters: general - 
reddening - metallicity 



1 INTRODUCTION 

The formation and evolution scenarios of the Milky Way Galaxy 
still remain among the most important unso lved problems in con- 
temporary astrophysics dPerrett et al7ll2002h . One promising way 
for us to better understand, and to possibly make progress towards 
addressing these issues is by studying globular clusters (GCs). GCs 
are generally considered the fossils of the galactic formation and 
evolution processes, since the y form during the ve ry early stages of 
their host galaxy's evolution dBarmbv et al. 1120001) . GCs are very 
dense, gravitationally bound spherical systems of several thousands 
to more than a million stars. They can be observed out to much 
greater distances than individual stars, so that they can be used, and 
are in fact the ideal tracers, to study the properties of extragalac- 
tic systems. The most distant GC systems that have been studied 
to dat e are those in the Coma cluster; for instance, Bau m et aT71 
dl995h presented GC counts in the br ight elliptical ga laxy NGC 
4881 (at a distance of ~ 108 Mpc: iBaum et al. Ill995l) based on 
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Hubble Space Te lescope (ffST)AVFPC2 observations. Using data 
of similar quality , Kavelaar s et al. Harris et al. 1( 12000) and 

IWoodworth & H arris (2000) published a series of papers on the 
GCs in NGC 4874 and IC 4051, the central cD galaxy and a giant 
elliptical galaxy in the Coma cluster, respectively. 

Located at a d istance of about 780 kpc dStanek & Garnavich I 
ll998j :[Macri 200l]), M31 is the nearest and largest spiral galaxy in 
the Local Group of galaxies. It contain s 337 confirmed GCs and 
688 GC candidates (Gallet i"et al. II2004T) . thus allowing us access 
to a much larger number of GCs than in our own Galaxy. How- 
ever, despite the difference in GC nu mbers, from the ob servational 
evidence collected so far (see, e.g.. iRich et al. I [2005b . the M31 
GCs and their Galactic counterparts reveal some striking similari- 
ties jFusi Pecci et al.lll994l;lDiorgpvski et al. ll997l : lB~armbv et al.1 



I2002X More recently. iKim et al. I d2007l) embarked on a new sys- 
tematic wide-field CCD survey of M31 GCs, and found 1 164 GCs 
and GC candidates - of which 559 are previously known GCs and 
605 newly-found GC candidates; based on survey data from the 
Canada France Hawaii telescope and Wide Field Camera on the 
Isaac Newton telescope. lHuxor 1 12007) combined his detailed anal- 
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ysis of the M3 1 GC system with recent results based on the galaxy's 
stellar halo, and concluded that M3 1 and the Milky Way are rather 
more similar than previously thought. Therefore, studying the prop- 
erties of the GCs in M31 may not only improve our understand- 
ing of the formation and structure of our nearest large neighbour 
galaxy, but also that of our own Galaxy. However, there are also sig- 
nificant difference s0 between the GC systems of the Milky Way and 
M31: in particular, M3 1 has a much larger population of GCs than 
the Milky Way (see lKim et al. |2007| . and references therein); there 
are populations of "faint fuzzies" and extended GCs in the outer 
halo of M31 that are not se en in the Milky Way dHuxor et al. 120051 ; 
iMackev et al.ll2006l, l2007h . and there is a population of young to 
intermediate-age G Cs in M3 1 that are (again) not seen in th e Milky 
Way again (see, eg jBeaslev et al. ll2004l : |Puzia et al. 1120051) . 

A reliable estimate of the reddening, caused by dust contam- 
ination, is important for the study of the stellar population of a 
given GC, in order to obtain its intrinsic spectral energy distribu- 
tion. In general, on galaxy-wide scales dust tends to be distributed 
close the the galactic plane in galactic discs. Therefore, the disc 
clusters are most affected by extinction due to dust in the galac- 
tic disc (in addition to the effects due to Galactic foreground ex- 
tinction). A reliable estimate of the extinction caused by Galac- 
tic material along a given line of sigh t can b e obt ained from the 
redden ing maps of iBurstein & Heiles I d 19821) and ISchlegel et al.1 
(1998). However, a reliable estimate of the internal extinction in a 
given galaxy is not easy to obtain. More specifically, although there 
were some ways in which to deal with the problem of determining 
the reddening toward the GCs in M31 prior to the publication of 
iBarmbv et al. d2000l). some may n ot have been fully satisfactory. 
For example, Ivan den Berghldl969l) assumed a uniform reddening 
for the clusters when he pioneered integrat ed-light spectroscopy of 
M31 GCs, whereas live & Richterl dl985t) assumed GCs to have 
a single intrinsic colou r when using them as reddening probes ; 
some authors (see e.g. iKron & Mayalll Il960l: IVetesnikl 1 1 962a ; 
lHarrislll974l ; Ivan den Berghlll975l; baiaia & Gergelvlll977l) as- 
sumed that the halo GCs in M31 were only affected by Galactic 
foreground extinction, based on which they then averaged the in- 
trinsic colours of the halo GCs and subtracted this value from the 
observed colours to obtain the reddening values for all the GCs (see 
details in IBarmbv et alE oOO). There are also two other ways to de- 
termine the reddening of M31 clusters, which seem r easonable but 
they h ave only been applied to a handful of clusters: iFrogel et al. I 
(1980) estimated the individual reddening for 35 clusters using 
the reddening-free para meter Qk, based on unp ublished spectro- 
scopic data by L. Searle uCrampton et al. Idl985l) used the intrinsic 
colours of the 35 clusters obtained by Frog el et al. I Jl98Ch to cali- 
brate (B — V)o as a function of spectroscopic slope parameter S of 
the continuum between ~ 4000 and 5000 A, and then determined 
the i ntrinsic colours for abo ut 40 GCs and GC candidates. 

IBarmbv et al. I d2000t) presented a new catalogue of pho- 
tometric and spectroscopic data of M31 GCs, and determined 
the reddening for each individual cluster using correlations be- 
tween optical and infrared colours and metallicity, and by defin- 
ing various "r e ddenin g-free" parameters based on this catalogue. 
IBarmbv et alTI d2000h found that the M31 and Galactic GC ex- 
tinction laws (see their table 6), and the M31 and Galactic GC 
colour-metallicity relations are similar to each other. They then 



1 Huxor (2007) suggests that the primary difference between the Galaxy 
and M31, and between their GC systems in particular, is likely due to the 
more vigorous recent merger history of M3 1 . 



estimated the reddening to M31 objects with spectroscopic data 
using the relation between intrinsic optical colours and metallic- 
ity as determined for Galactic clusters. For objects without spec- 
troscopic data, they used the relationships between the reddening- 
free para meters and certain intr insic colours, based on the Galactic 
GC data. IBarmbv et al~l d2000l) compared their results with those 
in the literature and confirmed that their estimated reddening val- 
ues are reasonable, and quantitatively consistent with previous de- 
terminations for GC s across the entire M31 disc. In particular, 
IBarmbv et afl feoOO) showed that the distribution of reddening val- 
ues as a function of position appears reasonable in that the objects 
with the smallest reddening are spread across the disc and halo, 
while the objects with the largest reddening are concentrated in the 
galactic disc. The re ddening values for M31 cl usters obtained by 
Barm by et al. Id200d) are widely used (see, e.g., Jjiang et al. 1120031 : 
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To study the metal abundance properties of GCs can help 
us understand the formation and enrichment processes of their 
host galaxy. For example, if galaxies form as a consequence of 
a monolithic, dissipative and rapid collapse of a single massive, 
nearly-spherical spinning gas cloud in which the enrichment time- 
scale is shorter than the collapse time, the halo stars and GCs 
should show large-s cale metallicity gr adients l Eggen et al. 19621 : 
IBarmbv et al. |2000|) . On the other hand. lSearle & Zinn Id 1978b pre- 
sented a chaotic scheme for the early evolution of a galaxy, in which 
loosely bound pre-enriched fragments merge with the main body 
of the proto-galaxy over a significant period, in which case there 
should be a more homogeneous metallicity distribution. At present, 
most galaxies are thought to have formed through a combination of 
both of these scenarios (see also Section 4.5) 

HST provides a unique tool for studying GCs in ex- 
ternal g alaxies. For example, ba s ed on data from th e HST 
arc hive. iGebhardt & Kissler-Patig I dl999h . Larsen et al.l d200lh 
and Kundu & Whitmore I feooih showed that many large galaxies 
possess two or more subpopulations of GCs that have quite differ- 
ent chemical co mpositions (see also IWest et al. 1 120041) . Recently, 
IPeng etal.ld200rl presented the colour distributions of GC systems 
for 100 early-type galaxies from the ACS Virgo Cluster Survey, and 
found that, on average, galaxies at all luminosities appear to have 
bimodal or asymmetric GC colour/metallicity distributions. The 
presence of colour bimodality among these old GCs indicates that 
there have been at least two major st ar-forming mechanisms in the 
(early) histories of massive ga laxies dWest et aL~ll2004l ; |Peng et aTl 
120061; iBrodie & Straderll2006l) . 

Based on the newest photometric and spectral data, in this pa- 
per we determine reliable reddening values for 443 GCs and GC 
candidates (the largest GC sample in M31 used to date), and we 
also determine the metallicities for 209 GCs and GC candidates 
without spectroscopic observations. We then perform a statistical 
analysis using this GC sample. We describe the photometric and 
spectroscopic data for the M31 GCs in Section[2] In Section[3] we 
determine and analyse the reddening of our M31 GC sample. Sec- 
tion [4] is devoted to our statistical analysis. Finally, the discussion 
and conclusions of this paper are presented in Section|5] 



2 DATABASE 
2.1 The sample 

iGalleti et al. I d2004l) present the final Revised Bologna Catalogue 
of M31 GCs including 337 confirmed GCs and 688 GC candi- 
dates which compose our primary sample. From a comparison with 
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iBarmbv et all feOOd) and lPerrett et al. I d2002l) . 89 candidates turn 
out not to be GC s, and these w ere thus removed from the sam- 
ple. More recently. [Huxor 1 d2007h provided a further revision of the 
Bologna Catalogue, including a number of additional clusters in the 
halo of M31. Because he only provides photometry for these new 
clusters in two filters, we will not include these in our final sample, 
however (see below regarding the photometric requirements of the 
method we adopt in this paper). 

2.2 The optical and near-infrared photometric data 

The source o f the photome t ric dat a utilised in this paper is from the 
catalogue of Gall eti et al. ] J2004I) . i.e. the updated Bologna Cata- 
logue with homogenised optical (UBVRI) photometry collected 
from the most recent phot ometric references available in the lit- 
erature. In this c atalogue, iGalleti et all d2004l) used the UBVRI 
photometry from lBarmbv et al. I ( 120000 as a reference to obtain the 
master catalogue of pho tometric measureme nts as homogeneously 
as possible. In addition. I Galleti et al. I J2004) identified 693 known 
and candidate GCs in M31 using the 2MASS database, and in- 
cluded their 2MASS JHK a magnitudes, tr a nsformed to the CIT 
photometric system felias et al. II 1982L 1 1983h . IGalleti et al.l d2004l) 
compile a final table including the UBVRIJHK S photometric 
data for the 337 confirmed and 688 candidate GCs in M31 (their 
table 2), which is the photometric database we use in this paper. 



2.3 Spectroscopic metallicities 

iHuchra et all Jl99lh obtained spectroscopy of 150 M31 GCs with 
the Multiple Mirror Telescope (MMT). The system they used has a 
resolution of 8-9 A and enhanced blue sensitivity. Their observa- 
tions extend to the atmospheric cut-off at 3200 A, since many of the 
strongest and most metallici ty-sensitive spectral features of interest 
are in the UV (see details in lBrodie & Huchrall99fj|) . lHuchra et al. I 
determined the metallicities for these 150 GCs using six 
absorption-line indices f rom integrated cluster spectra, following 
IBrodie & Huchra I dl990l) 

IBarmbv et al. I yOOOJ) determined the metallicities of 61 M31 
GCs and GC candidates using the Keck Low Resolution Imag- 
ing Spectrometer (LRIS) and the MMT Blue Channel spec- 
trograph. They computed th e abso r ption-line indices usi ng the 
method of IBrodie & Huchra I Jl990b . IBarmbv et al~l d200d) com- 
bined the measured ind i ces ba sed on the metallicity calibration 
from [Brodie & Huchra ] dl99d) , in order to determine metallici- 
ties. Their metallici ties were found to be consist ent with those of 
IHuchra et all dl99ll) . Finallv. frjlirmbv et al. I d200d) compiled a cat- 
alog ue of spectroscopic metallicities for 188 M31 GCsq 

IPerrett et alTI d20020 determined the metallicities of 201 GCs 
in M31 using the Wide-Field Fibre Optic Spectrograph (WYF- 
FOS) at the 4.2m William Herschel Telescope. They calculated 12 
absorp tion-line indices following the method of IBrodie & Huchra I 
(1990). By comparison of the line indi ces with the previ - 
ously published M31 G C [F e/H] values of i Bonoli etal.l dl987l), 
IBrodie & Huchra I Jl990t) . and lBarmbv et al. I d2000l) , |Perrett et all 
d2002l) found that the line indices of the CH (G band), Mg b and 
Fe53 lines best represent the metal abundances of their observed 
targets. They determined the final metallicities of these GCs based 
on an unweighted mean of the [Fe/H] values calculated from the G 
band, Mg b, and Fe53 line strengths. 
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There are some GCs in our sample for which the metal- 
licities were determined in two or three of the studies men- 
tioned above. To use all the met allicities as coheren tly as pos- 
sible, we use th e metallicities of | Perrett et al. I J2002T ) whenever 
available, since IPerrett et al. 1 d2002n~ present the largest (homoge- 
neous) sam ple of M31 GC meta llicit ies. For the metallicit ies de- 
termined by IHuchra etaTl dl99ll) and lBarmbv et al. I d200d) . there 
is on ly one GC in comm on, object B328. IHuchra et al. f dl99ll) 
and lBarmbv et al. 1 d2000h determined its metallicity at [Fe/H] = 
-1.22±0.80 and [ Fe/H] = -1.51±0.2 8, respectively. We adopt 
the metallicity from Barmb v" ;t al. Ijioool) , given its smaller associ- 
ated uncertainty. 

Overall, we obtained metallicities for 295 M31 GCs, which 
we list in Table Q] We will refer to these data as our spectroscopic 
metallicity catalogue (hereafter SMCat). We will use the SMCat to 
perf orm our statistical anal ysis. 

IBarmbv et al. ] d2000h and IPerrett et al. ] d2002l) determined 
the GC metallicit i es usi ng the metallicity calibration defined in 
Brodi e & Huchra 1 dl99d) . Therefore, all three metallicity de termi- 
nations are on the same [Fe/H] system. IPerrett et al. I d2002L their 
fig. 7) show convincingly that there are no systematic offsets among 
these three sets of metallicity determinations. 



3 REDDENING DETERMINATIONS 

As already discussed in the introduction, there are several ways of 
dealing with t he problem of determ ining the reddening towards the 
M31 clusters. [Barmbv et al. I J2000h determined the largest number 
of reliable reddening values for M31 GCs using correlations be- 
tween optical and infrared colours and metallicity, and by defining 
various "reddening-free " parameters based on their catalogue of 
multicolour photometry. Barm bv et al. 1 d200d) compared their re- 
sults with those in the literature and confirmed that their estimated 
reddening values are reasonable, and quantitatively consistent with 
previous determinations for GCs across the entire M31 disc. Be- 
low, we will d etermine more M31 G C reddening values based on 
the method of Barmby et al. ( 2000), and on the Revised Bologna 
Catalogue dGalleti et al. II2004I) . which is the newest and the most 
comprehensive multicolour catalogue available at present. 

3.1 Constructing the correlations based on Galactic clusters 

In this section, we will construct the correlations between opti- 
cal colours and metallicity by defining various "reddening-free" 
par ameters (hence f orth r eferred to as Q parameters), follow- 
ing Barmb v" ;t al. I (120 00). based on the infrared photometry of 
IBrodie & Huchra I d 19900 and on the updated Galactic GC cata- 
logue of lHarris 1 1 19961. up dated in February 2003; hereafter H03). 
Brodie & Huchra (1990) presented infrared photometry for 23 
Galactic GCs. H03 lists the reddening values, metallicities, and op- 
tical colours of 150 Galactic GCs. 

First, we performed linear regressions of intrinsic optical 
colours versus metallicit y. We use the Galactic GCs with E(B — 
V) < 0.5 mag, following lBarmbv et al. I d2000h : 



(X - Y) = a [Fe/H] + b 



E(B - V) = 



E(B-V) 
E(X - Y) 



[(X -Y)-(X- Y) ] 



(1) 
(2) 



where (X — Y) represents any colour, and (X — Y)q represents 
the relevant intrinsic colour obtained based on the reddening val- 
ues listed in H03. The reddening ratio can be determined from the 
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Table 1. Spectroscopic metallicities of the M31 GCs collected in this paper. 



Name [Fc/H] source Name [Fc/H] source Name [Fe/H] source Name [Fe/H] source 
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z 
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Table 1. Continued. 



Name 


[Fe/H] 


source 


Name 


[Fc/H] source 


Name 


[Fc/H] 


source 


Name 


[Fe/H] source 


B204 


-0.80 ±0.17 


1 


B207 


-0.81 ±0.59 1 


B208 


-0.84 ±0.04 


1 


B209 


-1.37 ±0.13 1 


B210 


-1.90 ±0.32 


1 


B213 


-1.02 ±0.11 1 


B214 


-1.00 ±0.61 


1 


DAO065 


-1.80 ±0.36 1 


DAO066 


-1.82 ±0.26 
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-1.87 ±0.39 1 


B217 


-0.93 ±0.14 


1 
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-1.19 ±0.07 1 
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-0.01 ±0.57 
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-1.21 ±0.09 1 
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-1.29 ±0.04 
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-0.93 ±0.95 1 
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-1.80 ±0.05 1 
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B372 


-1.42 ±0.17 1 
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-1.90 ±0.67 
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-1.86 ±0.66 1 
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-1.79 ±0.72 


1 


B483 


-2.96 ±0.35 1 


B484 


-1.95 ±0.59 


1 


B378 


-1.64 ±0.26 1 


B380 


-2.31 ±0.45 


1 


B382 


-1.52 ±0.27 1 


B386 


-1.62 ±0.14 


1 


B289D 


-1.71 ±0.63 1 


B292D 


-0.47 ±0.54 


1 


G327 


-1.88 ±0.06 1 


B391 


-0.55 ±0.59 


1 


B400 


-2.01 ±0.21 1 


BA11 


-1.14 ±0.61 


1 







NOTE - The spectroscopic metallicities used in this paper are from Per rett et al. I j2002l) (source=l). lBarmbv et al. I j2000j) (source=2) and lHuchra et al. I 

Jl99ll) (source=3). 



Galactic extinction law o f 1CardeilietaTl ll989). The fit results with 
correlation coefficients r > 0.8 are listed in Table [2] We us e bi- 
sector linear fits, as described by I Akritas"& Bershady l l 1996b . be- 
cause we are not only interested in the case where metallicity is 
used to predict colour, but also in the re verse case where co lour is 
used to predict metallicity (see details in lBarmbv et al. 1 2000). 

Next, we construct relationships between Q parameters and 
intrinsic colours, to estimate the reddening for clusters without 
spectroscopic metallicities. The Q parameters are defined as 



Q X r^(X-Y)-§^n(Y-Z) 

= {X- Y) - |^ ~ g (Y - Z) ,0) 

where X, Y and Z refer to photometric magnitudes in any filter. 
The relation between an intrinsic colour and the Q parameter is 

(X - Y) = aQxYZ + b , (4) 

(X - Z) = aQxYZ + b , (5) 

and 

(Y - Z) = aQxYZ + b (6) 

The fit results with correlation coefficients r > 0.8 are listed 
in Tabled 

Fig. 1 shows the relationships of a few representative fits be- 
tween Q-parameter and colour for Galactic GCs, randomly selected 
from the set of relations included in Table 3. 



3.2 Tests of the methods using heavily reddened Galactic 
clusters 

In this section, we will test the methods adopted to derive the 
reddening values, using the heavily reddened Galactic GCs with 
E(B — V) > 0.5 mag from H03 (which were not used to construct 
the calibrations discussed above). 

Based on Eqs. (l)-(6) and the correlation parameters from Ta- 
bles|2]and[3] we can determine the reddening values for these highly 
reddened Galactic GCs. 

For each of the two methods we averaged all values of E(B — 
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Figure 2. Comparison of the reddening values between H03 and this paper, 
for Galactic GCs. The error bars represent the uncertainties on the redden- 
ing values derived in this paper only; H03 does not provide uncertainty 
estimates. 



V) to produce one final value of E(B — V) per method. The stan- 
dard deviation of the average value of E(B — V) is taken as its 
error for each method. The result of the comparison is shown in 
Fig-IU from which we can see that the results are encouraging. The 
average offset between E(B — V) from the Q-parameter method 
and the H03 value is 0.01 ± 0.01 mag; for the colour-metallicity 
method, the average offset is 0.00 ± 0.01 mag. It is clear that the 
two data sets agree very well. 



3.3 Reddening values of the M31 clusters 

iBarmbv et ail < |2000|) showed that the M3 1 and the Milky Way red- 
dening laws are the same within the observational errors. There- 
fore, in this section, we will determine the reddening values for the 
M3 1 clusters and cluster candidates based on the calibrated colour- 
metallicity (C-M) and Q-parameter relations for the Milky Way 
from Tables [2] and [3] The metallicities are fro m the SMCat and th e 
optical and infrared photometric data are from lGalleti et al. I d2004f) . 
as discussed in Sections 2.2 and 2.3. For each object, we average 
all reddening values obtained using the various C-M and Q-colour 
relations, to get one value for the reddening. The uncertainty in the 
reddening value thus derived is calculated as the standard deviation 
of the resulting reddening values. 

We determined the reddening values for all M31 GCs and GC 
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Table 2. Colour-metallicity relations for Galactic GCs 









(X - Y) = 


o[Fe/H] + b 
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Table 3. Q-parameter-colour relations for Galactic GCs 
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candidates with sufficient data, a total of 658 objects. However, 
some reddening values are not reliable, such as those based on only 
one C-M or Q-colour relation, and those w ith large reddening er - 
rors. In order to maintain consistency with Bar mbv et al. I (2000), 
we adopted the rules followed by these authors who rejected red- 
dening values a E{B -v) / E(B - V) > 0.5 for E(B - V) > 
0.15 mag and a E (B-v) I E( B-V) > 1.0 for E (B - V) < 
0.15 mag. However, following iBarmbv et al~l (2000), we empha- 
size that the rules adopted here for rejecting reddening values are 
quite arbitrary. In total, 443 reliable reddening values were deter- 
mined in this paper, which are listed in Table|4] Columns 1, 3, 5, 
7, and 9 list the names of the GCs, using the nomenclature adopted 
bv lGalletietaI1 ( l2004h . 

From Sections 3.1 and 3.2 we find that the reddening values 
for the Galactic GCs obtained from different relations are internally 
consistent. However, for M31 GCs this is not always the case. For 
some GCs and GC candidates the reddening values, based on dif- 
ferent relations, are inconsistent. Reasons for this may include that 
for Galactic GCs, the photometric data are accurate, but for some 
M31 GCs and GC candidates (particularly the fainter objects) this 
may not be the case. Therefore, when we average the reddening 



values for each object, we reject reddening values that are clearly 
statistical outliers: these are defined as those values that differ from 
the mean value for a given object by more than la. 

Fig. [3] shows the distribution of the reliable reddening values 
listed in Table|4] From Fig.|3]we find that slightly more than half of 
the reddening values are E(B — V) < 0.2 mag. The distribution 
of the 443 reliable reddening values has a mean of E(B — V) = 
0.281 

a 14 w ^ n a stan dard deviation of a — 0.17 mag, compar ed 
with E(B - V) = 0.22; o = 0.19 mag of lBarmbv et ail bOOOl) . 

Fig. [4] shows the reddening values as a function of posi- 
tion. The l arge ellipse rep resents the boundary of the M31 disc 
defined by iRacine I ( 1199 11) and the small ellipses on the north- 
western and southeastern sides of the major axis are the D25 
diameters of the M31 companion galaxies NGC 205 and M32, 
respectively. The distribution appears reasonable in that the ob- 
jects with low reddening values are spread across the disc and 
halo, while those with high reddening values are mainly concen- 
trated in the galactic disc. However, from Fig. [4] we can also see 
that a substantial number of objects outside the "halo" boundary 
have E(B — V) > 0.1 mag, i.e., greater than the Galactic fore- 
ground reddening in the direction of M31, as estimated by many 
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4. Reliable reddening values of GCs and GC candidates in M3 1 . 



Name 


E(B - V) 


Name 


E(B - V) 


Name 


E(B - V) 


Name 


E(B - V) 


Name 


E(B - V) 


G001 


1 4- 09 
U. -LU ZC U . UZ 


G002 


A f)K _L A 01 
U.UU zz U . U -L 


B290 




BA21 


u.oo zn u.ui 


B412 


9fi 4- 09 
u.zu zc. u.uz 


B413 


A AQ -U fl 04 


B 134D 


n 4^ 4- n oq 

u.tu zn u.uy 


B291 


n -4- n n9 
u.uo zn u.uz 


B138D 


n 9Q 4- n (14 

U.ZO ^ U.U^L 


B 1 40D 


4^ 4-011 
\j.*±i) in u. 1 1 


B141D 


u.io nz u.uy 


B142D 


A r:o _L a oi 
U.UO zz U.UO 


B144D 


u.oo ^ u.u^ 


B147D 


n qq 4_ n ncr 
u.oo ^ u.uo 


B295 


1 4- 01 

U. .LU ZZ U. U .L 


B 148D 


n 04 -r- n 01 
u.u*± in u.uo 


B 149D 


A ftQ _|_ A on 

u.uy zn u.uu 


B 150D 


u.oz ^ u.uz 


B416 


(1 ^4 4- (1 (1^ 

U.O^L ^ U.UO 


B 152D 


x.O 4- 1 Q 
u.uz m u. iy 


B418 


9fi 4- OQ 

u.zu m u.uy 


B 154D 


a OQ 4- 07 
u.uy it u.u i 


B 156D 


n 4^ -4- n no 
u.4o ^ u.uy 


B420 


n in 4 (i (ii 
u.ou ^ u.uo 


B 157D 


1 f; 4 n 01 
u. iu in u.ui 


B422 


A 1 A _L A OS 
U. -LU ZZ U.UO 


B 162D 


A c;q _L A 07 
U.UO IT U.UI 


B 163D 


U. -LU U.UO 


R9Q8 


n 1 4- n (19 

U. -LU □! U.UZ 


B424 


f^Q 4- 01 
u.uy in u.uo 


B165D 


A Ol _|_ A CY7 

U.Z _L ZH U.UI 


B166D 


9fi 4- 04 

U.ZU zz u.ut 


B301 


(1 17+ (1(19 
u.li zn u. u^ 


B167D 


n 9i 4- n ni 

u.^o ^ u.uo 


B427 


e;i 4- o 01 
U.UO ZZ u.uo 


B169D 


a n7 4- fl 01 
U.U 1 zn u . uo 


B170D 


91 4- 1 1 
U. ZO ZZ U. -L J_ 


B302 


n in 4- n (11 

u. lu zn U. U L 


B428 


(1 ^4 4- (1 (1^ 
u.o*± ^ u.uo 


B172D 


11 -1-0 0^ 

U. -L _L ZZ U.UU 


B430 


a i a _L a nc; 
u.iu zn u . uu 


B173D 


14 4- 01 

u.ot m U.UO 


B175D 


(1 (IS 4- (1 (I 11 ! 
u.uo ^ u.uo 


B303 


n 1 4 4- n ofi 

U. L<± ^ u.uu 


B177D 


OS 4- 09 
u.uo ZZ u.uz 


B304 


n 07 + n ni 

U.UI _l_ U.UI 


B433 


o (Si +00^ 

U.UI _1_ U.UU 


B306 


n 42 + (1 02 

U.rlii _l_ U.Ui 


B435 


1) fi7 + fl ITi 

U.UI _l_ u.uu 


B307 


OS + 02 

U.UO _1_ U.Uii 


B178D 


a i A 4- n 09 

U. -Lrfc Z1Z U.UZ 


B309 


17 4- 04 
U. -L I zz u.ut 


B310 


OQ 4- 01 
u.uy ± u. u .L 


B181D 


Ifi 4- OQ 
u.ou zz u.uy 


B31 1 


9Q 4- 09 
u.zy zz u.uz 


B438 


A QO _L A 07 


B312 


1 fi 4- 01 

U. -LU X U.UI 


B183D 


Ifi 4- 0^ 
u.ou m u.uu 


B314 


OS 4- 0^ 
u.uo m u.uu 


B313 


91 4- o 09 
u.z_l m u.uz 


B315 


a n7 4- 09 

U.U 1 ZZ U . UZ 


B001 


4- 09 
U.ZU zz u.uz 


B316 


91 4- o 01 

U.ZrJ. U.UO 


B317 


Oil 4- o 09 

U . -L J. ZZ U . UZ 


B186D 


11 4- o O 1 ^ 
U.OO ZZ u.uu 


B440 


A QO _|_ A 1 K 
U.OZ ZH U. -LO 


B003 


1 Q 4- 09 
u. -i_y zn u.uz 


B188D 


c;fi 4- Ofi 
u.uu m u.uu 


B 190D 


9fi 4- 01 
u.zu m u.uo 


B004 


07 4- 09 
u.ui in u.uz 


B005 


A OQ _L A 00 

U.ZO IH u . uz 


B325 


1 4 4- 09 
u. it zz u . uz 


B328 


1 o 4- 09 

U. -LU ± U. UZ 


B192D 


A^ 4- 09 
U.IU ZZ u.uz 


B330 


n 4- o 01 

U. O -L ZZ U.UO 


R004H 


A c;7 _L A 07 
U.UI ZH U.U i 


R006 


OQ 4- 09 
u.uy zn u.uz 


B 194D 


c;i 4- 04 
u.uo m u.ut 


B447 


14 4- O 1 1 
u.ot m u. -LO 


B244 


97 4- 01 
U.Z I ZZ u.uo 


ROOQ 


A 1 Q _L A 09 

U. -LO in u.uz 


roi n 

DU1U 


99 4- o 01 
u.zz zn U.UI 


B01 1 


11 4- o 01 
u. 1 1 m u.ui 


B012 


19 4-001 
u. -Lz m u.ui 


B196D 


1 Q 4- Ofi 

u. iy in u.uu 


B245 


1 Q7 _L A 07 
l.OI IH U.UI 


B448 


O 1 ^ 4- 01 
u.uu zn U . U -L 


B013 


1 1 4- o 09 

U. -LO ± U. UZ 


B014 


ifi 4- o 09 
U.OU ZZ u.uz 


B197D 


iq 4- o f)^ 
u.oy zz u.uu 


B335 


A «CC _L A 09 
U.UU ZIZ U.UZ 


B449 


i 97 4- OS 

J..ZI ZIZ U.UO 


B015 


c;o 4- 09 
u.uu ± u.uz 


B016 


10 4- 09 
u.ou ± u.uz 


B450 


94 4- 1 

U.Z'l ZC U. -LU 


B337 


Ofi 4- 09 
U.UU ZZ u . uz 


BO 17 


97 4- 09 
u. z i in u . uz 


B018 


90 4- 01 

U.ZrU ± U . U -L 


B019 


90 4- 01 

U.ZU ± U. U -L 


B020 


19 4-001 

U. -LZ ZZ U.U -L 


B338 


a i a 4- fl 09 

U. -L'i ZIZ U.UZ 


B021 


9fi 4- 09 
u. zu in u . uz 


B022 


04 4- 01 
U.UI ± u.uo 


B339 


1 fi 4- 01 

U. -LU ZZ U.UO 


B023 


19 4- 01 

U.OZ ZZ U. U -L 


B453 


10 4- 09 
U.OU ZH U.UZ 


B024 


01 4- 09 
u.uo m u.uz 


V031 


11 4- 09 
u.oo m u.uz 


B025 


90 4- 01 
u.zu m u.ui 


B202D 


IS 4- 07 
u.oo m u.u i 


B027 


A Ol ±0 01 
U.Z.L ZIZ U.UI 


B026 


1 k 4- o 09 
u. .LU in u . uz 


B028 


99 4- 09 


B020D 


99 4- Ofi 

U.ZZ ZIZ u.uu 


3029 


1 9 4- 01 

U. .LZ ZZ U. U .L 


DUjU 


ZlS -r- 01 

U.4:0 ZH U.UO 


D\)J 1 


11 4- 09 
u.oo m u.uz 


DUJZ 


49 4- 09 
u.tz m u.uz 


R4S6 


O 19 4- 04 
u.oz m u.ut 


B203D 


Ifi 4- 04 
u.ou m u.U'i 


B033 


1 4 4- 09 
U. -L^i ZZ U . UZ 


B034 


1 Q 4- o 01 

u . J_y zn u . u .L 


B457 


1 4 4- 09 

U. Lt ^ U . UZ 


B036 


1 z. _l o 09 

U. -LU ZZ U . UZ 


B204D 


AS 4- 99 


B025D 


A KQ _L A 01 

u . uo in u . uo 


B037 


1 Ol _L o 01 

_l .z _l zn u . uo 


B038 


97 4- 01 

U.Zr 1 U. U -L 


B039 


IS 4- 09 
u.oo ± u.uz 


B205D 


7S 4- Ofi 
U. 1 o ZZ u . uu 


B041 


07 4- 01 

U.UI zz U.UO 


B042 


fil 4- 01 

U.U .L zn U . U .L 


B044 


11 4- 01 

U.OO ± U. U -L 


B343 


Ofi 4- 01 

U.UU ZZ U. U -L 


B045 


1 S 4- 01 

U. LO ZZ U. U .L 


B046 


1 Q 4- 01 

u. iy zn u.uo 


B207D 


11 4- OS 
u.oo m u.uo 


B048 


1 Q 4- 09 
u. xy m u.uz 


L3UH / 


OQ 4- 09 
u.uy m u.uz 


R04Q 


1 fi 4- 09 
u. iu m u.uz 


B050 


a OA 4- n 01 

U. Z'i ZH U. U -L 


B051 


14 4- o 09 
u.o'-t m u.uz 


B052 


91 4- o 04 
u.zo zn u.ui 


B054 


91 4- o 09 
U.ZO ZZ u.uz 


B056 


17 4-0 01 

U.L 1 ZZ U.U-L 


B057 


a OQ 4- 09 
u.uy m u . uz 


B058 


1 1 4- 01 

U . -LO m U. U .L 


B059 


9Q 4- 01 
u. zy ^ u. u .L 


B060 


91 4- o 09 
U.Z J. ZZ u.uz 


B061 


14 4- 09 
U.OI ZZ u.uz 


B062 


9fi 4- 01 

U. ZU m u . uo 


B063 


40 4- 01 
u.iu m u. u .L 


B064 


174-OOI 

U. -L 1 ZZ U. U .L 


B065 


10 4-001 

U . -LU m U . U .L 


3344 


11 4- o 09 

U. .L _L ZZ U.UZ 




94 4- 01 
U.Z^ zn U.UO 


RftfiS 

DUOO 


IS 4- 01 
u.oo m u.uo 


B257 


i 17 4-0 01 
i . -L i m u.uo 


B461 


c;s 4- 07 
u. uo m u.u / 


B070 


19 4-0 04 
u. iz m u.U'i 


B073 


All ±001 
U. -L -L ZH U . U -L 


B072 


11 4- o Ofi 
u.oo m u.uu 


B074 


1 Q 4- o 01 
u. xy zn u . u ± 


B075 


10 4- 04 
u.ou ZZ u.ut 


3077 


Q7 4- o 01 
u. y i zz u.uo 


B078 


n 44. 4- 07 
u. it m U.UI 


B079 


c;e 4- o 01 
u.uo m u.uo 


B081 


11 4- o 09 

U. -L .L ZZ U.UZ 


B345 


1 4- 09 

U. -LU ZZ U . UZ 


B462 


1Q 4- 04 
u.oy zz u .\) L ± 


B082 


fi9 4- 01 
U.UZ m U.UO 


B083 


19 4 n (19 
u . lz m u. uz 


B084 


9fi 4- 04 

U.ZrU ZZ u.ut 


B085 


1 4 4- 09 

U. Lt ZZ U.UZ 


3086 


1 c; 4- o 01 

U. .L U ZZ U. U .L 


B346 


71 4. a oi 
U. 1 J. m u . uo 


B088 


4fi 4- 01 
u.iu m u . u .L 


B092 


1 9 4- 09 

U. -LZ ^ U. UZ 


B347 


1 4 4- 09 

U. Lt ^ U.UZ 


3348 


9 c; 4- 04 

U.ZU ZZ U.U'i 


B093 


10 4" 09 

U.UU _1_ U.Ui 


B094 


07 + 02 

U.UI _l_ U.Ui 


B095 


41 4- o 04 


B096 


9fi 4- 02 

U.ZiU _l_ U.Ui 


3098 


OS + 02 

U.UO _1_ U.Uii 


B463 


a ii 4- 07 
u. oo zz U.UI 


B097 


9Q 4- 01 
u. zy ^ u . u .L 


B099 


1 fi 4- 01 
U . -LU ZZ U.UO 


B350 


1 4- 09 

U. -LU ZZ U.UZ 


3100 


AS 4- OS 
U.^iO zz u.uo 


B101 


17 + 002 

U.l 1 _l_ U.Ui 


NB46 


fi9 4- 07 

U. _l_ U.UI 


B103 


1 Q 4- o 02 

U.-LJ7 _l_ U.Ui 


NB70 


QQ 4- 01 

U.Uf _l_ u.uu 


B464 


27 + 01 

U.iil _1_ u.uu 


B105 


U. -Lrt _l_ U.UI 


B106 




B107 


n 28 + 02 

U.iiO _l_ U.Ui 


B109 


n (18 + 02 

U.UO _l_ U.Uii 


Bill 


n (18 4- 02 

U.UO _1_ U.Uii 


B110 


90 4- 09 

U.ZU zz u.uz 


B260 


fi7 4- 09 
U.UI ± u . uz 


B112 


1 4 4- 09 

U. Lt ^ U . UZ 


B117 


04 4- 01 

U.UI ± U . U -L 


Bl 15 


19 4-0 01 
U. .LZ ZZ U . U .L 


Bl 16 


fi9 4- 09 

u.uz m u . uz 


NB64 


4fi 4- 1 ^ 

U.IU ± U. -LU 


B118 


99 4- Ofi 
U.ZiZ ^ u.uu 


B119 


1 4 4- 09 

U. Lt ^ U.UZ 


3351 


1 c; 4- 09 

U. .L U ZZ U.UZ 


B352 


a i a 4- n 09 
U. J-'i zn u.uz 


NB25 


fifi 4- o 1 7 

U.UU ± U. -L 1 


B122 


QQ 4- 09 
U.OU ZZ u.uz 


B123 


10 4- 01 
U.OU ZZ u.uo 


3125 


O^t 4- 09 
u.uu zz u.uz 


B127 


0Q 4- 09 
u.uy m u . uz 


B354 


O^i 4- 09 
u.uu ± u.uz 


B128 


1 e; 4- o 01 

U. -LU ZZ U. U .L 


B129 


i 1 fi 4- Ofi 

L.LU ZZ U.UU 


NB39 


AS 4- 09 
U.'iO zz u.uz 


B130 


Ifi 4- 01 
U.UU _l_ U.UI 


B131 


1 9 4- 04 

U. -LZi _l_ U.Urt 


B134 


01 + 02 

U.UU _l_ U.Ui 


B135 


07 4- o 01 

U.Zi 1 _l_ U.UI 


3355 


07 + Ofi 

U.Ul _1_ u.uu 


B136 


Ifi 4- 04 

u.ou zz u.ut 


B137 


40 4- 09 
U.IU ± u.uz 


B217D 


1 a 4- o 01 

U . -LO ZZ U . U .L 


B141 


19 4- 01 

U.OZ ^ U . U -L 


3143 


4- 09 
u.uu zz u.uz 


B144 


0^ -r- 09 

u.uo in u.uz 


B219D 


O 49 4- O 09 
u.tz m u.uz 


B146 


Ofi 4- 04 
u.uu m u.ut 


B266 


QS 4- OQ 
u.yo in u.uy 


3148 


1 7 4- 09 
u. 1 1 m u.uz 


R790D 


a 07 4- 04 
U.UI zz u.ut 


B149 


izL 4- 09 
u.ot ± u.uz 


B221D 


c;i 4- o OS 

U. UO ZZ u.uo 


B467 


97 4- o 09 

U.Z 1 ZZ u.uz 


3150 


9S 4- 07 

U.ZO ZZ U.UI 


B223D 


0.20 ±0.05 


B151 


0.32 ± 0.01 


B152 


0.18 ± 0.01 


B356 


0.31 ±0.01 


B153 


0.05 ±0.01 


B154 


0.15 ±0.03 


B468 


0.27 ±0.03 


B357 


0.12 ± 0.02 


B155 


0.20 ±0.02 


B156 


0.10 ±0.02 


B158 


0.14 ±0.00 


B159 


0.36 ± 0.04 


B226D 


0.63 ± 0.00 


B161 


0.17 ±0.01 


B162 


0.25 ±0.03 


B163 


0.14 ±0.01 


B358 


0.06 ± 0.01 


B164 


0.12 ± 0.03 


B165 


0.12 ±0.01 


B228D 


0.18 ±0.02 


B167 


0.03 ±0.02 


B168 


0.54 ±0.05 


B169 


0.59 ± 0.04 


B170 


0.10 ±0.02 


B272 


0.57 ±0.04 


B171 


0.11 ±0.01 


B172 


0.18 ± 0.02 


DAO062 


1.11 ± 0.17 


B173 


0.40 ± 0.04 


B174 


0.32 ±0.01 


B177 


0.18 ±0.03 


B176 


0.10 ± 0.01 


B178 


0.12 ± 0.01 


B179 


0.10 ±0.01 


B180 


0.14 ±0.01 


B181 


0.23 ±0.01 


B231D 


0.07 ±0.02 


B182 


0.25 ± 0.01 


B185 


0.11 ±0.01 


B184 


0.21 ±0.03 


B470 


0.40 ±0.08 


B187 


0.35 ± 0.02 


B471 


0.62 ± 0.04 


B189 


0.04 ±0.03 


B190 


0.10 ±0.02 


B192 


0.31 ±0.02 


B194 


0.07 ±0.02 


B193 


0.11 ± 0.01 


B472 


0.13 ±0.00 


B195 


0.12 ±0.00 


B196 


0.26 ± 0.04 


B235D 


0.49 ± 0.03 


B197 


0.19 ± 0.02 


B199 


0.10 ±0.02 


B201 


0.04 ±0.02 


B202 


0.26 ±0.02 


B203 


0.16 ± 0.02 


B204 


0.12 ± 0.01 


B361 


0.11 ±0.01 


B237D 


0.32 ±0.05 


B205 


0.14 ±0.01 


B206 


0.13 ± 0.01 


B238D 


0.40 ±0.07 


B208 


0.13 ±0.02 


G260 


0.30 ±0.05 
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(u-B) 





E(B-V) 

Figure 3. Distribution of the reddening values of the M31 GCs and GC 
candidates obtained in this paper. 



authors (see, e.g., van den Bergh 1969; McClur e & Racine |[l969l ; 
iFrogel et al. 1119801) . In fact, iBarmbv et al. I j2OO0t> also noted this 
phenomenon. They suggested a number of plausible explanations, 
which include that (i) this could be caused by the large uncertain- 
ties inherent to the method; or (ii) the assumption that the M31 
halo clusters are subject to only foreground reddening is somewhat 
dubious. 

IVetesnikU1962bh ana lysed the photometry of M31 GCs pub- 
lished bv lVetesnikl dl962al) . He assumed that the halo clusters are 
only affected by foreground Galactic extinction. He derived a mean 
intrinsic colour of (B — V)o = 0.83 mag from the objects in the 
halo of M3 1 and calculated the colour excess for each object. Sub- 
sequently, he studied the reddening distribution of M31 GCs on 
either side of the major axis and found that the GCs on the north- 
western side of the disc are either intrinsically redder, or suffer from 
more significant extinction, than those on the southeastern side of 
the disc. 

Based on three homogeneo us and independent p hotometric 
data sets for M3 1 GC candidates, live & Richter I d 1985b examined 
the differential reddening towards these objects, and found that the 
GCs on the northwestern side of the disc are redder than those on 
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Table 4. Continued. 



Name 


E(B - V) 


Name 


E(B - V) 


Name 


E(B - V) 


Name 


E(B - V) 


Name 


E(B - V) 


B239D 


0.38 ± 0.12 


B209 


0.10 ± 0.01 


B211 


0.10 ± 0.01 


B212 


0.13 ± 0.01 


B213 


0.15 ± 0.02 


B214 


0.05 ± 0.02 


B215 


0.21 ± 0.04 


B362 


0.65 ± 0.20 


G268 


0.26 ± 0.01 


B217 


0.12 ± 0.01 


G270 


0.67 ± 0.20 


B218 


0.14 ± 0.01 


B219 


0.05 ± 0.03 


B243D 


0.03 ± 0.02 


B220 


0.05 ± 0.02 


B245D 


0.52 ± 0.03 


B221 


0.23 ± 0.02 


B222 


0.05 ± 0.02 


B224 


0.13 ± 0.02 


B473 


0.20 ± 0.04 


B225 


0.10 ± 0.01 


B226 


1.08 ± 0.06 


B247D 


0.51 ± 0.08 


B227 


0.37 ± 0.04 


B228 


0.13 ± 0.01 


B229 


0.07 ± 0.02 


B230 


0.15 ± 0.01 


B365 


0.19 ± 0.02 


B231 


0.15 ± 0.03 


B232 


0.14 ± 0.01 


B233 


0.17 ± 0.01 


B281 


0.12 ± 0.02 


B250D 


0.54 ± 0.06 


B252D 


0.26 ± 0.04 


B234 


0.11 ± 0.02 


B366 


0.05 ± 0.02 


B474 


0.53 ± 0.05 


B475 


0.16 ± 0.03 


B235 


0.11 ± 0.01 


B256D 


0.69 ± 0.06 


B476 


0.08 ± 0.05 


B236 


0.07 ± 0.05 


B258D 


1.20 ± 0.09 


B237 


0.14 ± 0.02 


B260D 


0.20 ± 0.03 


B478 


1.00 ± 0.13 


B370 


0.34 ± 0.01 


B238 


0.11 ± 0.02 


B239 


0.09 ± 0.01 


B261D 


0.27 ± 0.06 


B263D 


0.25 ± 0.07 


B240 


0.13 ± 0.00 


B286 


0.67 ± 0.02 


B479 


0.64 ± 0.11 


B266D 


0.75 ± 0.15 


B372 


0.20 ± 0.02 


B373 


0.10 ± 0.01 


B375 


0.29 ± 0.03 


B481 


0.52 ± 0.08 


B377 


0.16 ± 0.02 


B270D 


0.25 ± 0.02 


B483 


0.08 ± 0.06 


B378 


0.14 ± 0.02 


B379 


0.15 ± 0.01 


B273D 


0.29 ± 0.04 


B274D 


0.23 ± 0.04 


B380 


0.06 ± 0.02 


B381 


0.17 ± 0.02 


B275D 


0.18 ± 0.03 


B486 


0.17 ± 0.02 


B277D 


0.34 ± 0.08 


B382 


0.10 ± 0.02 


B278D 


0.42 ± 0.04 


B384 


0.04 ± 0.02 


B385 


0.06 ± 0.05 


B386 


0.21 ±0.01 


B283D 


0.16 ±0.04 


B288D 


0.58 ±0.06 


B387 


0.12 ±0.02 


B489 


0.17 ±0.04 


B289D 


0.23 ±0.05 


B490 


0.30 ±0.02 


G325 


0.12 ±0.09 


B389 


0.27 ±0.08 


B293D 


0.27 ±0.06 


G327 


0.18 ±0.01 


B295D 


0.62 ±0.12 


B296D 


0.06 ± 0.04 


B297D 


0.30 ±0.07 


B298D 


0.70 ±0.11 


B299D 


0.39 ±0.00 


B300D 


0.74 ±0.15 


B393 


0.14 ±0.02 


B492 


0.26 ±0.03 


B302D 


0.25 ± 0.02 


B304D 


0.67 ±0.04 


B493 


0.09 ±0.07 


B494 


0.57 ±0.07 


B307D 


0.45 ±0.08 


B308D 


0.12 ± 0.07 


B495 


0.34 ±0.08 


B396 


0.09 ±0.01 


B310D 


0.34 ±0.03 


B313D 


0.11 ±0.04 


B314D 


0.12 ± 0.03 


B317D 


0.69 ±0.03 


B319D 


0.44 ±0.06 


B320D 


0.50 ±0.06 


B324D 


0.11 ±0.06 


B398 


0.16 ± 0.03 


B399 


0.03 ±0.02 


B400 


0.21 ±0.02 


B326D 


0.31 ±0.12 


B328D 


0.65 ±0.09 


B329D 


0.16 ± 0.03 


B330D 


0.70 ±0.02 


B331D 


0.50 ±0.06 


B332D 


0.33 ±0.13 


B402 


0.16 ±0.03 


BA11 


0.06 ± 0.03 


B334D 


0.42 ±0.05 


B338D 


0.55 ±0.09 


B339D 


0.61 ±0.05 


B403 


0.07 ±0.02 


B340D 


0.23 ± 0.06 


B405 


0.14 ±0.02 


B508 


0.10 ±0.08 


B343D 


0.18 ±0.03 


B344D 


0.07 ±0.04 


B345D 


0.08 ± 0.05 


B346D 


0.19 ±0.06 


B407 


0.16 ±0.02 


B347D 


0.16 ±0.05 


B348D 


0.42 ±0.04 


B349D 


0.21 ± 0.01 


NB63 


0.90 ±0.05 


NB16 


0.41 ±0.08 


NB50 


0.93 ±0.21 











the southeastern side. Therefore, live &RichterH 1985b concluded 
that the northwestern side of the disc of M3 1 is nearer to us than its 
southeastern side. 

With the large sample of M31 clusters (and candidates) at 
hand, we can now also examine the reddening distribution of M3 1 
objects on either side of the major axis, and calculate the mean 
E(B - V). 

Following IPerrettetal. I fc002h and lHuchra et al. I (1991), we 
use the X, Y plane to indicate the position of the GCs. The X co- 
ordinate is the position along the major axis of M31, where positive 
X is in the northeastern direction, while the Y coordinate is along 
the minor axis of the M31 disc, increasing towards the northwest. 
The relative coordinates of the M3 1 clusters are derived by assum- 
ing standard geometric parameters for M31. We adopted a central 
position for M31 at a = 00 h 42 n M4 3 .30 and dp = ±41°16'09". 
(J2000.0) following lHuchra et al. I ( fl99lh and lPerrett et al. I j2002h . 
Formally, 



X = A sin (9 + 73 cos 
Y = -Acos8 + Bsint 



(7) 
(8) 



where A — sin(a — an)cos<5 and B = sin5cos<5o — cos(a — 
Qo) cos 8 sin So. W e adopt a po sition angle of 6 — 38° for the 
major axis of M31 dKentlll989h . Fig. [5] shows graphically the de- 
pendence of the average reddening on the distance from the major 
axis. The error bars represent the standard deviations of the means. 
It is clear that the reddening on the northwestern (NW) side of 
the disc is much greater than that on the southeastern (SE) side. 
The mean reddening on the northwestern and southeastern sides is 
E(B - V) = 0.33 ± 0.02 and 0.24 ± 0.02 mag, respectively. 




Y (arcmin) 

Figure 5. Reddening distribution dependence on the distance from the ma- 
jor axis of M3 1 . The dashed lines indicate the mean reddening values for 
GCs on the two sides. 



Below, we will check our resulting reddening values by study- 
ing the distribution of the colours of the M31 clusters and cluster 
candidates as a function of the projected distance, Y, from the ma- 
jor axis. If our reddening values are correct, the distribution of the 
intrinsic colours of the M31 clusters should be symmetric, even if 
the distribution of the observed colours is asymmetric. The left- 
hand panel of Fig. [6] shows the distribution of the mean colours of 
the GCs and GC candidates binned in 5.5 arcmin intervals in Y. 
The error bars represent the standard deviations of the means. It 
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Figure 4. Map of the confirmed and candidate GCs in M31. The large ellipse is the M31 disc/halo boundary as defined by Racine] J 199 lb ; the two small 
ellipses are the D25 isophotes of NGC 205 (northwest) and M32 (southeast). 



is clear that the colour distribution is asymmetric. The right-hand 
panel of Fig.[6]shows the distribution of the mean intrinsic colours 
of the GCs and GC candidates binned in 5.5 arcmin intervals in Y. 
Again, the error bars represent the standard deviations of the means. 
Obviously, the distribution of the mean intrinsic colours is nearly 
symmetric. On the southeastern side, the mean intrinsic colour is 
(B — V)a= 0.64±0.03, while on the northwestern side the mean 
intrinsic colour is (B — V)o=0.66±0.02. 



4 M31 GLOBULAR CLUSTER METALLICITIES 

It is evident that the metallicity distribution of a galaxy's GC sys- 
tem can provide important clues as to the process and conditions 
relevant to galaxy formation. Previous studies of the M31 GC 
system have presented some important information. For example, 
sig ns of bimodali t y were found in the M3 1 metalli c ity distribution 
by IHuchra et al. I ll99ll) . lAshman & Birdl Jl993l). iBarmbv et al I 
( 2000) and lPerrett et al. I d2002h . In addition. IHuchra et al. I dl99lh 
found that the metal-rich clusters in M3 1 appear to form a central 
rotating disc system. Based on the current largest sample including 
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Y (arcmin) 



Y (arcmin) 



Figure 6. (B — V) colour and intrinsic (B — V)o colour versus distance 
from the major axis of M31 for M31 GCs and GC candidates. 



321 velocities. IPerrett et al. l J200l) performed a more comprehen- 
sive investigation into the kinematics of the M3 1 cluster system and 
showed that the metal-rich GCs appear to constitute a distinct kine- 
matic subsystem that demonstrates a centrally concentrated spatial 
distribution with a high rotation amplitude, but it does not appear 
to be significantly flattened. This is consistent with a bulge popula- 
tion. In this section, we will examine the distribution of GC metal- 
licities in M3 1 using the largest number of GCs and GC candidates 
available to date. We will include the metallicities determined based 
on GC colours. 

4.1 Colour-derived metallicities 

In total, 23 1 of the GCs and GC candidates in our sample of 443 ob- 
jects with reliable reddening values have no spectroscopic metallic- 
ities. Therefore, we will determine their metallicities based on the 
C-M relation for these objects. From the colou r exces ses of GCs 
with spectroscopic metallicities, Bar mbv et all (2000) found that 
the M3 1 and Galactic GC C-M relations are consistent. Therefore, 
in this section, we will determine metallicities for 23 1 GCs and GC 
candidates from their intrinsic colours by applying the Galactic C - 
M relation. We use bi-sector linear fits ( Akritas & Bershadv 8 19961) 
to determine the metallicity as a function of colour, and average 
the resulting metallicities over the available colours for each object 
(except for some significantly deviating points that are most likely 
due to inaccurate photometric data, i.e. those that differ from the 
mean value by more than la, as justified above). As in the red- 
dening determination, the standard deviations of the metallicities 
from individual colours are used as the error estimates. There are 
209 GCs (and candidates) with metallicity determinations, which 
are listed in Table [5] 

4.2 Comparison of spectroscopic and colour-derived 
metallicities 

In this section, we will test the process outlined in Section 4. 1 by 
applying it to the clusters with spectroscopic metallicities; this in- 
cludes all GCs for which the metallicities can be determined based 
on the C-M relation fits. The results of our comparison are shown 
in Fig. [7] the mean metallicity offset (spectroscopic minus colour- 
derived metallicity) is 0.039 ± 0.022 dex. From Fig.|7]we can see 
that there is no evidence of a bias in the prediction of the metallici- 
ties. The offsets for all clusters are less than 0.7 dex, and the offsets 
for 4 clusters are greater than 0.5 dex. Since some of these metal- 
licity differences are substantial, we have carefully double checked 



our data and results. In fact, the large spread in m etallicity space is 
also evident from fig. 10 of Barmbv et al. (2000). However, upon 
close examination of the data, we found that there are 4 objects 
that should not be included in the calculation of the offset between 
spectroscopic and colour-derived metallicities, i.e. objects B068, 
B075, B159 and B219. For B068, we determined 14 colour-derived 
metallicities, and these colour-derived metallicities have a small 
scatter. Their mean value is [Fe/H] — —0.81 ± 0.03 dex. How- 
ever, the spectral metallicity is [Fe/H] = —0.29 ± 0.59 dex; we 
suspect that the spectral metallicity determination may be incor- 
rect. For B075, there are only three colour-derived metallicities, the 
mean value of which is [Fe/H] = —1.71 ± 0.08 dex. On the other 
hand, the spectral metallicity is [Fe/H] = —1.03 ± 0.33 dex, so 
that we think that more photometric data is needed to determine 
the colour-derived metallicities more accurately. For B159, there 
are also only three colour-derived metallicities. Their mean value 
is [Fe/H] = -2.33 ±0.11 dex, whilst the GC'sspectral metallicity 
is [Fe/H] = —1.58 ± 0.41 dex. Therefore, we also think that more 
photometric data is needed to determine its colour-derived metal- 
licity more accurately. Finally, for B219 we determined 14 colour- 
derived metallicities, which exhibit a small scatter. The mean value 
of these metallicities is [Fe/H] = —0.62 ± 0.05 dex. However, 
the spectral metallicity is [Fe/H] = —0.01 ± 0.57 dex, while we 
also suspect that this spectral metallicity may be problematic. Ex- 
cept for these four clusters, the mean metallicity offset (spectro- 
scopic minus colour-derived metallicit y) is 0.028±0.022 dex , com- 
pared with 0.020 ± 0.021 found by (Barmbv et al. I |2000) based 
on a smaller GC sample. This bias in the metallicity determina- 
tion may come from large errors in either the colour- or spectro- 
scopically determined metallicities, or both. In addition, the corre- 
lations between optical colours and metallicity, which are used to 
determine the colour-derived metallicity for M31 clusters, are con- 
structed based on the Galactic GCs, which may have introduced a 
small (but likely insignificant) bias. In the following analysis, we 
have substracted this offset from all of our colour-derived metal- 
licities. Fig. [8] shows the metallicity distributions of the spectro- 
scopic and colour-derived samples (cf. Fig. [7}. Here, we now use 
a Kolmogorov-Smirnov (KS) test to demonstrate whether the two 
distributions in Fig. [8] are the same. We determined a value of 
C mal = 0.061 for these two samples (which do not include the 
four objects noted above); D mal is defined as the maximum value 
of the absolute difference between two cumulative distribution 
functions. The probability of obtaining a value of D mal = 0.061 
is 80.8%. It is clear that the KS test supports the similarity of the 
metallicity distributions of the spectroscopic and colour-derived 
samples. 

4.3 Metallicity distribution 

The metallicity distribution of the M31 clusters has been in - 

199lh, 



Huchra et al. 



and Perrett et all 



vestigated in pr e vious s tudies, including _ 
Ashm an & Bird I dl993l) . iBarmbv et al~l bOOO 
(2002). With the current largest GC and GC candidate sample at 
hand, we will now reanalyse the M31 GC metallicity distribution. 
Including the metallicities determined based on the C-M fits, our 
sample includes a total of 504 metallicities. 

Fig. [9] shows the metallicity distribution of the M31 GCs dis- 
cussed in this paper, along with a similar figure for the Milky Way 
GCs (from H03), for comparison. To assess whether the combi- 
nation of spectral and colour-derived metallicities has an effect 
on the results, we consider three data sets in our analysis of the 
metallicity distribution of the M31 objects: Set 1 contains all ob- 
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Table 5. New metallicity estimates for M31 GCs and GC candidates without spectroscopic observations. 



Name 


[Fe/H] 


Name 


[Fe/H] 


Name 


[Fe/H] 




Name 


[Fe/H] 


Name 




[Fe/H] 




B290 


— 1.07 ± 0.07 


BA21 


—2.51 ± 0.11 


B412 


—0.80 ± 0. 


.02 


B413 


— 1.57 ± 0.08 


B134D 


—2 


.41 ± 0. 


.23 


B291 


— 1.12 ± 0.04 


B138D 


—0.36 ± 0.04 


B140D 


— 1.57 ± 0. 


.12 


B141D 


— 1.07 ± 0.08 


B142D 


—2 


.59 ± 0. 


.24 


B144D 


— 1.62 ± 0.09 


B148D 


—1.93 ± 0.21 


B149D 


—2.21 ± 0. 


.25 


B150D 


—2.52 ± 0.07 


B416 


— 1 


.34 ± 0. 


.07 


B152D 


—2.55 ± 0.09 


B418 


— 1.19 ± 0.10 


B154D 


—0.53 ± 0. 


.83 


B156D 


—2.58 ± 0.14 


B420 


— 


.63 ± 0. 


.07 


B157D 


—0.09 ± 0.08 


B422 


— 1.97 ± 0.18 


B162D 


—2.53 ± 0. 


.17 


B163D 


0.20 ± 0.06 


B424 


—2 


.09 ± 0. 


.10 


B165D 


— 1.05 ± 0.13 


B166D 


— 1.02 ± 0.08 


B167D 


—2.34 ± 0. 


.08 


B427 


— 1.57 ± 0.08 


B169D 


— 


.10 ± 0. 


.04 


B428 


—2.04 ± 0.09 


B172D 


—2.51 ± 0.12 


B175D 


0.31 ± 0. 


.09 


B177D 


— 1.32 ± 0.01 


B433 


—2 


.45 ± 0. 


.14 


B435 


— 1.76 ± 0.12 


B178D 


—0.90 ± 0.24 


B309 


—2.03 ± 0. 


.26 


B181D 


—2.21 ± 0.20 


B438 


—2 


.42 ± 0. 


.31 


B186D 


—2.08 ± 0.14 


B440 


—0.41 ± 0.09 


B003 


—2.08 ± 0. 


.07 


B188D 


— 1.84 ± 0.14 


B325 


— 1 


.77 ± 0. 


.08 


B192D 


—0.46 ± 0.03 


B330 


— 1.98 ± 0.06 


B004D 


— 1.17 ± 0. 


.05 


B194D 


— 1.97 ± 0.12 


B447 


— 1 


.73 ± 0. 


.12 


B244 


— 1.50 ± 0.21 


BOM 


—0.55 ± 0.13 


B197D 


—0.69 ± 0. 


.05 


B450 


—0.89 ± 0.05 


B022 


— 1 


.64 ± 0. 


.07 


B339 


—0.90 ± 0.04 


B202D 


—0.74 ± 0.04 


B020D 


—0.76 ± 0. 


.08 


B032 


— 1.09 ± 0.03 


B456 


— 1 


.71 ± 0. 


.29 


B203D 


—2.03 ± 0.08 


B457 


1 st\ I r\ o "i 

— 1.60 ± 0.21 


B204D 


— 1.18 ± 0. 


.10 


B025D 


— 1.28 ± 0.10 


B205D 


—2 


.35 ± 0. 


.21 


B207D 


—0.27 ± 0.08 


B052 


a 1 n I n 1 ^7 

0.12 ± 0.17 


B060 


— 1.87 ± 0. 


.06 


B062 


—0.47 ± 0.11 


B067 


—2 


.55 ± 0. 


.04 


B257 


—2.05 ± 0.82 


B461 


—2.56 ± 0.07 


B070 


— 1.66 ± 0. 


.10 


B078 


—0.56 ± 0.15 


B079 


— 


.85 ± 0. 


.03 


B345 


— 1.40 ± 0.06 


B462 


—2.28 ± 0.34 


B084 


—0.76 ± 0. 


.07 


B346 


— 1.70 ± 0.07 


B347 


— 1 


.71 ± 0. 


.03 


B348 


— 1.38 ± 0.07 


B463 


— 1.46 ± 0.18 


B099 


— 1.03 ± 0. 


.06 


B 1 00 


—2.21 ± 0.10 


B101 


— 1 


.17 ± 0. 


.02 


NB46 


— 1.48 ± 0.03 


NB70 


—2.48 ± 0.04 


B464 


—0.44 ± 0. 


.09 


Bill 


— 1.50 ± 0.03 


B260 


— 


.36 ± 0. 


.10 


NB64 


—2.12 ± 0.46 


B 1 18 


— 1.64 ± 0.10 


B351 


— 1.60 ± 0. 


.05 


NB25 


—0.31 ± 0.08 


B123 


— 1 


.58 ± 0. 


.01 


B128 


—0.92 ± 0.05 


NB50 


—2.23 ± 0.19 


B136 


—2.39 ± 0. 


.08 


B217D 


—2.36 ± 0.04 


B266 


—2 


.80 ± 0. 


.15 


B220D 


—2.69 ± 0.13 


B221D 


— 1.56 ± 0.15 


B150 


—0.76 ± 0. 


.08 


B223D 


—0.23 ± 0.08 


B468 


—2 


.16 ± 0. 


.12 


B155 


r\ O A 1 A AO 

—0.84 ± 0.03 


B226D 


—2.01 ± 0.19 


B162 


—0.70 ± 0. 


.05 


B228D 


0.27 ± 0.14 


B168 


— 


.12 ± 0. 


.21 


B169 


—2.56 ± 0.06 


B172 


—0.87 ± 0.03 


DAO062 


—2.13 ± 0. 


.14 


B173 


— 1.86 ± 0.47 


B177 


— 


.88 ± 0. 


.10 


B181 


— 1.10 ± 0.03 


B231D 


—0.12 ± 0.07 


B470 


—2.18 ± 0. 


.23 


B187 


— 1.72 ± 0.04 


B471 


—2 


.18 ± 0. 


.09 


B189 


0.18 ± 0.11 


B194 


—1.56 ± 0.05 


B195 


— 1.48 ± 0. 


.63 


B196 


—1.94 ± 0.08 


B202 


— 1 


.84 ± 0. 


.11 


B361 


— 1.61 ± 0.02 


B237D 


—0.78 ± 0.26 


G260 


—2.45 ± 0. 


.06 


B239D 


— 1.67 ± 0.13 


B215 


— 1 


.21 ± 0. 


.03 


G268 


— 1.36 ± 0.01 


B243D 


— 1.28 ± 0.07 


B245D 


—2.88 ± 0. 


.09 


B473 


—2.17 ± 0.16 


B247D 


— 1 


.90 ± 0. 


.23 


B227 


— 1.28 ± 0.08 


B250D 


—0.98 ± 0.14 


B252D 


—2.83 ± 0. 


.09 


B474 


—2.12 ± 0.10 


B256D 


—2 


.37 ± 0. 


.13 


B476 


—0.03 ± 0.13 


B236 


— 1.01 ± 0.17 


B258D 


—2.47 ± 0. 


.08 


B260D 


— 1.46 ± 0.15 


B478 


—2 


.69 ± 0. 


.01 


B261D 


-2.45 ± 0.19 


B263D 


-0.85 ± 0.07 


B286 


-1.67 ±0. 


.11 


B479 


-0.36 ± 0.13 


B266D 


—2 


.34 ± 0. 


.10 


B481 


-1.45 ±0.13 


B270D 


-2.28 ±0.19 


B273D 


-1.01 ±0. 


.03 


B274D 


-0.08 ±0.11 


B275D 


-1 


.81 ±0. 


.13 


B277D 


-0.83 ±0.06 


B278D 


-2.46 ±0.19 


B385 


-0.86 ± 0. 


.14 


B283D 


-1.55 ±0.17 


B288D 


-2 


.58 ±0. 


.23 


B489 


-0.04 ±0.10 


B490 


0.08 ± 0.07 


G325 


0.12 ±0. 


.12 


B389 


-0.35 ±0.08 


B293D 


-2 


.57 ±0. 


.11 


B295D 


-2.22 ±0.13 


B296D 


-0.91 ±0.16 


B297D 


0.10 ±0. 


.08 


B298D 


-2.35 ±0.11 


B299D 


-1 


.91 ±0. 


.10 


B300D 


-2.52 ±0.11 


B393 


-1.41 ±0.05 


B492 


-1.06 ±0. 


.27 


B302D 


-1.30 ±0.05 


B304D 


-2 


.41 ±0. 


.12 


B493 


1.07 ±0.28 


B494 


-1.54 ±0.05 


B307D 


-1.64 ±0. 


.10 


B308D 


-2.43 ±0.14 


B495 


-0 


.35 ±0. 


.05 


B396 


-1.87 ±0.09 


B310D 


-0.93 ±0.05 


B313D 


-1.00 ±0. 


.16 


B314D 


-0.18 ±0.06 


B317D 


-2 


.47 ±0. 


.13 


B319D 


-1.86 ±0.17 


B320D 


-2.62 ±0.15 


B324D 


-2.27 ±0. 


.00 


B398 


-0.72 ±0.08 


B399 


-1 


.69 ±0. 


.09 


B326D 


-0.63 ±0.08 


B328D 


-1.67 ±0.06 


B329D 


-0.24 ±0. 


.04 


B330D 


-0.55 ±0.06 


B331D 


-1 


.50 ±0. 


.08 


B332D 


-0.65 ±0.09 


B402 


-1.18 ±0.06 


B334D 


-0.75 ±0. 


.06 


B338D 


-1.86 ±0.06 


B339D 


-1 


.51 ±0. 


.04 


B340D 


0.19 ±0.29 


B508 


-2.61 ±0.07 


B343D 


-0.49 ± 0. 


.05 


B344D 


-1.40 ±0.03 


B345D 


-0 


.39 ±0. 


.12 


B346D 


-0.37 ±0.19 


B347D 


0.00 ±0.15 


B348D 


-1.16 ±0. 


.26 


B349D 


-0.76 ±0.26 











jects for which metallicities have been determined from spectro- 
scopic observations; Set 2 contains the metallicities without spec- 
troscopic observations, which have been determined based on the 
C-M relationships in this paper; and Set 3 contains all metallicities 
from Sets 1 and 2. The mean metallicities of Sets 1, 2, and 3 are 
[Fe/H] = -1.202 ± 0.036, -1.414 ± 0.057, and -1.290 ± 0.032 
dex, respectively, i.e. comparable to (although not in all cases the 
same as) the value of [Fe/H] = —1.21 ± 0.03 dex obtained by 
IPerrett et al. I d2002h . and also comparable to the value of [Fe/H] = 
- 1.298 ± 0.046 dex obtained for the Milky Way GCs (H03). How- 
ever, we point out that the mean metallicities of Sets 1 and 2 are 
different at the 3a level. We investigated this result carefully, and 
found that there are many more objects of which the metallici- 
ties determined based on the C-M relationships in this paper are 
lower than [Fe/H] = —2.5, compared to the number of GCs for 
which metallicities were determined from spectroscopic observa- 
tions. Generally speaking, the objects without spectroscopic ob- 



servations are fainter than those with spectroscopic observations. 
Therefore, we should keep in mind that the objects for which the 
metallicities have been determined based on the C-M relationships 
are fainter than those with spectroscopic observations. The pho- 
tometry is more accurate for bright than for faint objects, however. 
For bright clusters, i.e., those with spectroscopic observations, the 
metallicities determined based on the C-M relationships are con- 
sistent with those determined from sp ectroscopic observations (see 
Fig. 8). In fact, Bar mbv et all (2000) did not use the very high or 
low metallicity values, [Fe/H] > 0.5 or [Fe/H] < —2.5 dex, ob- 
tained based on the C-M relationships in their paper when investi- 
gating the metallicity distribution. If we exclude these very high or 
low metallicity values, the mean metallicities of Sets 1 and 2 are 
the same. At the same time, we emphasize that there are no reasons 
for not using these low or high metallicities when investigating the 
metallicity distribution for M31 GCs. The asymmetric appearance 
of the metallicity distribution in Fig. [9] suggests the possibility of 
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Figure 7. Comparison of spectroscopic and colour-derived metallicities for 
M3 1 GCs with spectroscopic data; <r s indicates the uncertainty in the spec- 
troscopic metallicity and <r c is based on the C-M relations. 



that the metallicity distributions of both data sets show strong bi- 
modality. Thus, the KMM tests suggest that we can consider these 
two distributions of metallicity bimodal at the ~ 100 per cent con- 
fidence level. We also investigated the KMM results based on three 
groups and on heteroscedastic two-group fits. The results of these 
exercises are listed in Tables 7 and 8, which show that three-group 
and heteroscedastic two-group fits are also statistically acceptable. 
We point out that KMM tests assume Gaussian distributions, which 
may or may not be realistic here. In the following analysis, we in- 
vestigate the metallicity distribution for all samples in this paper, 
for which the homescedastic two-group fits may be more appropri- 
ate (see Tables 6, 7 and 8). We also apply the Dip test to confirm 
whether or not the metallicity distributions of our M31 GC sam- 
ples are multimodal. The Dip statistic is the maximum difference 
between the empirical distribution function, and the unimodal dis- 
tributi on function that minimizes the maximum di fference (see de- 
tails in Hartigan_&Fiartigan 1985; Hartigan 1985). The Dip statis- 
tic measures the departure of the sample's unimodality. If its value 
approaches zero, the samples are taken from a unimodal distribu- 
tion; if the Dip measure is positive, the samples come from a mul- 
timodal distribution. The Dip values for Sets 1, 2 and 3 are 0.024, 
0.025 and 0.018 with significance values of 71.2%, 55.43% and 
66.9%, respectively, which shows that the Dip values support the 
multimodality found from the KMM tests. We note that the Dip 
value for the (bimodal) metallicity distribution of the Galactic GCs 
is 0.039 with significance value of 90.5%. 
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Figure 8. Comparison of spectroscopic and colour-derived metallicity dis- 
tributions. 



bimodality. We use the KM M algorithm dMcLachlan & Basfordl 
1 19881 : lAshman et al. Ifl994h to search for bimodality in the metal- 
licity distribution. The input of the KMM algorithm includes the 
individual data points, the Gaussian group membership to be fitted, 
and starting values for the estimated means and common variance. 
In fact, for two-group homescedastic fitting, the KMM algorithm 
is insensitive to these input values and converges to the "correct" 
two-group fit even when its starting points are far from the true 
means and variance (see details in Ashm arTet al. |[T994h . We as- 
sumed the same variances for both groups in the metallicity distri- 
bution; in this homoscedastic fitting the p-value of the hypothesis 
test returned by the KMM algorithm adequately indicates that a 
two-group fit is an improvement with respect to a one-group fit. Ta- 
ble[6]lists the parameters returned by the KMM algorithm. It is clear 



4.4 Spatial distribution 

In the previous section, we investigated the metallicity distribution 
of the M31 GCs based on the KMM test. The KMM test allows 
us to distinguish between the metal-poor and metal-rich subsam- 
ples, i.e., the KMM test gives [Fe/H]=-1.18 dex as the dividing line 
between the metal-poor and metal-rich GCs. However, there are 
about 54 objects that exhibit intermediate probabilities of member- 
ship in both groups (0.5 < prob. < 0.6). Since it is difficult to 
decide unequivocally which group these objects belong to, we sim- 
ply adopted the dividing line between metal-poor and metal-rich 
GCs from the KMM test. 

Figure [10] shows the projected spatial distributions of the 
metal-poor and metal-rich GCs in M31. Using Eqs. (7) and (8), 
we obtained the distances to our sample clusters from the centre 
of M31. From Fig. [TUl it is clear that the metal-rich GCs in M31 
ar e more centrally conc entr ated, consistent with what was found 
bv lHuchraetal.Nl99lh and lPerrett etal.l d2002h . The metal-poor 
GCs appear to occupy a more extended halo, although also with 
a general concentration following the outline of the M31 disc. The 
latter may have been caused by selection biases, i.e., from Fig. 4 we 
can see that there are more clusters observed along the major axis 
than the minor axis. Fig. QT] shows the histograms of the metal- 
poor and metal-rich populations. A notable shortage of metal-poor 
clusters in the innermost radial bins can clear l y be s een, which is 
also consistent with the results of lPerrett et al. I ( I2002h . In the Milky 
Way, the metal-rich GCs reveal significant rotation a nd have his- 
torically been asso ciated with the thick-dis c system dZinn I | l985 : 
Armandrofflll989h: howev er, other studies dFrenk & White II 19821 ; 
Minniti ll995l ; IC6te ll999l ; lForbes et al. 1200 lh have suggested that 
metal-rich GCs within ~ 5 kpc of the Galactic Centre are more ap- 
propri ately associated with the M ilky Way's bulge and/or bar. In 
M31. lElson & Walterbosl dl988l) showed that the metal-rich clus- 
ters constitute a more highly flattened syste m than the metal-poor 
ones, and appear to have disc-like kinematics; lrluchra et al. I dl99lh 
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Figure 9. Metallicity distributions and homescedastic bimodal KMM tests of M31 GCs and GC candidates, subdivided by uncertainty, and Galactic GCs. 



Table 6. Results from the KMM homescedastic bimodality tests for the metallicities of the GCs in M3 1 and the Milky Way. 



Data Set 



[Fe/H] 



CT [Fc/H] 

[Fc/H]i [Fc/H] 2 CT [Fc/H] 



1 


-1.202 ± 0.036 


0.617 


-1.558 


-0.676 


0.438 


182 


113 


0.032 


2 


-1.414 ±0.057 


0.828 


-1.948 


-0.619 


0.508 


124 


85 


0.000 


3 


-1.290 ±0.032 


0.719 


-1.740 


-0.722 


0.510 


284 


220 


0.004 


MW 


-1.298 ± 0.046 


0.564 


-1.620 


-0.608 


0.306 


101 


47 


0.000 



showed that the m etal-rich GCs are pref erentially located close to 
the galactic centre. iHuchra et al71 ( 1 199 lh also showed that the dis- 
tinction between the rotation of the metal-rich and metal-poor clus- 
ters is most apparent in the inner 2 kpc. Therefore, they concluded 
that the meta l-rich clusters i n M31 appear to form a central rotating 
disc system. IPerrett et al. I {2002) performed a more comprehen- 



sive investigation into the kinematics of the M31 cluster system. 
They showed that the metal-rich M31 GCs appear to constitute a 
distinct kinematic subsystem that demonstrates a centrally concen- 
trated spatial distribution with a high rotation amplitude, but that 
does no t appear significantly fl attened, consistent with a bulge pop- 
ulation. Schrod er et al. I d2002h performed a maximum-likelihood 
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Table 7. Results from the KMM heteroscedastic bimodality tests for the metallicities of the GCs in M3 1 . 



Data Set 


[Fe/H] i 


[Fe/H] 2 


°"l,[Fc/H] 


CT 2,[Fc/H] 


m 


ni 


P 


1 


-1.487 


-0.574 


0.466 


0.401 


207 


88 


0.145 


2 


-2.424 


-1.192 


0.180 


0.743 


44 


165 


0.000 


3 


-1.794 


-0.786 


0.488 


0.534 


254 


250 


0.022 



Table 8. Results from the KMM homescedastic trimodality tests for the metallicities of the GCs in M3 1 . 



Data Set 


[Fe/H]i 


[Fe/H] 2 


[Fe/H] 3 


CT [Fc/H] 


ni 


ni 


U3 


P 


1 


-1.733 


-1.100 


-0.462 


0.374 


120 


114 


61 


0.128 


2 


-2.141 


-1.159 


-0.266 


0.399 


92 


78 


39 


0.000 


3 


-1.985 


-1.206 


-0.434 


0.419 


162 


237 


105 


0.012 
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0.05 
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Figure 11. Radial distribution of the metal-rich and metal-poor GCs and 
GC candidates in M3 1 . 



kinem atic analysis of 166 M31 clusters taken from lBarmbv et al~l 
(2000) and found that the most significant difference between the 
rotation of the metal-rich and metal-poor clusters occurs at inter- 
media te projected galactocentric radii. Particularly, Schrod er et al. I 
presented a potential thick-disc population among M31's 
metal-rich GCs. 



tween metallicity and projected radius, but most of his clusters 
were locat ed inside a radius of 50 arcmin; however, some authors 
(see. e.g..lHuchra etal.|[l982l : ISharov|[l988l ; iHuchra et al.1ll99ll; 
IPerrett et al. 112002b showed that there is evidence for a weak but 
measurable metallicit y gradient as a function of projected radius. 
iBarmbv et all (2000) confirmed the latter result based on their 
large sample of spectral and colour-derived metallicities. 

In Fig. [12] we show the metallicity of the M31 GCs as a func- 
tion of galactocentric radius based on our large cluster sample. It 
is clear that the dominant feature of this diagram is the scatter in 
metallicity at any radius. However, it is also true that a trend of de- 
creasing metallicity with increasing galactocentric distance exists, 
for both the metal-poor and the entire population. The slopes of the 
metal-poor subsample and for the entire sample are — 0.006±0.001 
and —0.007 ± 0.002 dex arcmin -1 , respectively, while for metal- 
rich sample, it is 0.000 ± 0.001 dex arcmin -1 . The latter can cer- 
tainly be considered as no metallicity gradient. In order to show 
this, we display the mean metallicity binned in 10 arcmin inter- 
vals in galactocentric radius. We can see that, within ~ 90 ar- 
cmin, the mean metallicity decreases with galactocentric radius 
for both the metal-poor and for the entire population. The error 
bars represent the standard deviations of the m eans. These results 
are in good agreement with lPerrett et al. I d2002h and lHuchra et al. I 
dl99ll) . Therefore, we can conclude that simple smooth, pressure- 
supported collapse models of galaxies by themselves are unlikely 
to fit M31. 



4.5 Metallicity gradient 

The presence or absence of a radial trend in the metallicity of 
a GC sample is an import ant test of galaxy f ormation theories 
jBarmbv et al. II2OO0I) . In the Eg gen et al. I d 19621) galaxy formation 
scenario, the ha lo stars and GCs should show large-s cale metal- 
licity g radients (lEggenetal.lll962l ; IBarmbv et al~ll2000h ; however, 
in the lsearle & Zinnldl978h scenario the expected m etallicity dis- 
tributi on is more homogeneous. For the Milky Way, Armandroff 
( 1989) provided some evidence that metallicity gradients with both 
distance from the Galactic plane and distance from the Galac- 
tic Centre are present in the disc cluster system. For M3 1 , there 
are some inconsistent conclusions, e.g., van den Bergh ] d 19691) 
showed that there is little or no evidence for a correlation be- 



4.6 Metallicity versus intrinsic magnitude 

The (possible) correlation between cluster mass (or luminosity) and 
metallicity is important in GC formation theory. It is generally be- 
lieved that if self-enrichment is important in GCs, the most mas- 
sive clusters could retain their metal-enriched supernova ejecta, so 
that the metal abundance should increase with cluster mass; the 
opposite is true if cooling from metals determines th e tempera- 
ture in the cluster-forming clouds ( IBarmbv et al~| [2000). The pos- 
sible self-enrichment of GCs has been s tudied in detail in some 
aspects (see details in lStrader et al. 112006]). However , the m odel of 
GC self-enrichment developed by IParmentier et al. I d 19991) is par- 
ticularly interesting in this context. In this model, cold and dense 
clouds embedded in the hot proto-galactic medium are assumed 
to be the progenitors of galactic halo GCs. Based on this model, 
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Figure 10. Spatial distribution of the metal-rich and metal-poor GCs in M3 1 . 
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Figure 12. Metallicity distribution versus projected radius for the entire and 
the metal-poor populations of the M3 1 GCs and GC candidates. 



studies, the metal-r ich (redder) GCs d id not show a correspond- 
ing trend (see also iBekki et al. I [2007), Ba sed on high -resolution 
cosmological simulations including GCs, IBekki et a Li f2007h in- 
vestigated the formation processes and physical properties of GC 
systems in galaxies, and found that luminous metal-poor clusters 
would develop a correlation between luminosity and metallicity if 
they originated from the nuclei of low-mass galaxies a t high red- 
shift. In fact, in the simulations of IBekki et al. I (2007), the "sim- 
ulated blue tilts" come from the assumption that luminous metal- 
poor clusters originate from the stellar galactic nuclei of the more 
massive nucleated galaxies exhibiting a luminosity - metall icity re- 
lation. It is therefore evident that, in Bek ki et al. I ( 120071) ■ galax- 
ies which experienced more accretion/merging events of nucleated 
low-mass galaxies are more likely to show a blue tilt. 

Fig. [T3] shows the diagram of GC metallicity versus dered- 
dened apparent magnitude. It is clear that there is no obvious 
trend o f me tallicity with luminos ity similar to that in lHuchra et al. I 
jl99ll) and Bar mbv et all (2000). Least-squares fits show no evi- 
dence for a relationship between luminosity and metallicity in our 
sample clusters. 



Parmentier & Gilmore 1 feooih suggested that the most metal-rich 
proto-GCs are the least massive ones. 

The HST provides a unique tool to study GC s in external 
^alaxi e s. Recently, us i ng the AC S onboard the HST, Harris et al. I 
20061) . iMieske et al. I d2006h and IStrader et al. I J20061) found that 
in giant ellipticals - such as M87, NGC 4649 and NGC 7094 (al- 
though not in NGC 4472) - luminous blue GCs reveal a trend 
of having redder colours, such that more massive GCs are redder 
(mo re metal-rich). This trend is referred to as the "blue tilt" (see 
also iBrodie & St rader 2006]). This blue tilt has b e en interpreted 
as a re sult of self-enrichment l lStrader et al. I|200"6T) . IStrader et al. I 
j2006h specula tively sug gested that these GCs once possessed dark 
matter haloes. ISpitler et al. I J2006I) subsequently found that this 
blue tilt is also present in the Sombrero spiral galaxy (NGC 4594) 
and may extend to l ess luminous GCs w ith a somewhat shallowe r 
slo pe than derived bvlrlarris et al. I J20061) and lStrader et al. I J20061) . 
As lSpitler et al. I ( feOOq) pointed out, the Sombrero galaxy provides 
the first example of this trend in a spiral galaxy, and in a galaxy 
found in a low-density galaxy environment. However, in these ACS 



5 DISCUSSION AND CONCLUSIONS 

In this paper, we have (re-)determined the reddening values for 
443 clusters and cluster candidates in M31, as well as metallic- 
ities for 209 sample objects without spect roscopic observations . 
We have followed the methods described bv lBarmbvetal.l ( l200rj) . 
who found that the M3 1 and Galactic extinction laws are the same 
within the observational errors, and that the M31 and Galactic GC 
C-M relations are also consistent with each other. The sample of 
spectroscopic and photometric data used in this paper is the newest 
and largest to date. The spectroscopic data were obtained from the 
most recent references currently available and the photometric data 
are from the most comprehensive catalogue of M3 1 clusters avail- 
able at present, which includes 337 confirmed GCs and 688 GC 
candidates. Using the metallicities of the largest sample of clusters 
and cluster candidates at hand, we studied the properties of the M3 1 
clusters. Our main conclusions are summarised below: 

(i) The reddening distribution shows that slightly more than 50 
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Figure 13. Metallicity as a function of the reddening-corrected magnitude 
Vb for our sample GCs. 



per cent of the GCs suffer from a reddening of less than E(B — 
V) = 0.2 mag, and the mean value is E(B-V) = 0.28±^ mag. 
The spatial distribution of E(B — V) indicates that the reddening 
on the northwestern side of the M3 1 disc is greater than that on the 
southeastern side, which is consistent with the conclusion that the 
northwestern side in nearer to us. 

(ii) The metallicity distribution of the M31 GCs is bimodal with 
peaks at [Fe/H] « -1.7 and -0.7 dex. 

(iii) The diagram of metallicities as a function of radius from the 
M31 centre shows a metallicity gradient for the metal-poor GCs, 
but no such gradient for the metal-rich GCs. 

(iv) The metal-rich clusters appear to constitute a centrally con- 
centrated spatial distribution; however, the metal-poor clusters tend 
to be less spatially concentrated. 

(v) There is no correlation between luminosity and metallic- 
ity among our M31 sample clusters, which indicates that self- 
enrichment is indeed unimportant for cluster formation in M3 1 . 



We reiterate that in using the method o f lBarmbvet all 12000), 
there are two major unavoidable assumptions (acknowledged by 
these authors), i.e. that in the Milky Way and in M31 both the 
extinction law and the intrinsic colours of the GCs are the same. 
The latter assumption seems reasonable, since there is no evi- 
dence that GCs in different galaxies have different intrinsic colours. 
Regarding the former assumption, there is inconsistent evidence 
as to whether or not this is a valid assumption. For example, 
Walterbos & Kennicutt ] dl988l) found that the extinction law in 
M31 is very similar to that in the Milky Way, by analysing the 
two major dust lanes on the near side of M31; however, several 
studies have suggested that the reddening in M31 a ppears to be pe- 
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■B)/E(B- 
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V) 
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1.01 ±0. 



0.5 dMassevetal.lll995l) 



compared to 0.72 for the same ratio in the Milky Way. Based on 
a larg e sample of GCs wit h optical and near-infrared photometric 
data, iBarmbv et al. 1 12000) demonstrated that the U - and AT-band 
extinction curve of M3 1 is consistent with that of the Milky Way, 
with total-to- selective extinction coefficient Rv = 3.1. In fact, the 
former assumption is plausible because in the M31 disc the compo- 
sition and size distribution of the large normal grains which dom- 



inate the d ust mass may be s imilar to those in the Milky Way (see 
for details. lXu & Heloulll996l) . 

As an example, we will discuss in some detail the reddening 
value of the M31 GC B037 (a.k.a. 037-B327), which is known to 
be an extremely red o bject. There are a few references that dis - 
cuss this GC, incl udin glBarmby et aTl d2002bh.lMa et al. I l2006al) . 
iMa et al. I J2006ch and lcohen I J2006h . iKron & Mavall I dl960h first 
noticed an extremely red colour in photographic (P) and visual (V) 
bands for B037, and determined its absorption to be Av = 3.90 
mag. Bas ed on the photom e tric data for M31 star clusters in U, B, 
and V of I Vetesnik U 1962ah , Ivetesnik I dl962bl) studied the redden- 
ing values for these objects and found that B037 was the most 
highly reddened in his sample, with E(B ~ V) = 1.28 mag 
(A v = 4.10 mag). ICrampton et al. I ( fl98a) calibrated (B - V)o 
as a function of spectroscopic slope parameter 5* of the contin- 
uum between ~ 4000 and 5000 A, and then determined the intrin- 
sic colours for about 40 GCs and GC candidates, including B037. 
ICrampton et al. I J 1985t> presented a reddening value for B037 of 
E{B — V) — 1. 48 mag. Armed with a large database of multi- 
colour photometrv. lBarmbv et al. 1 d2000h determined the reddening 
value for each individual M3 1 GC, including B037, using the cor- 
relations between optical and infrared colours and metallicity based 
on various "reddening-free" parameters, and derived E(B — V) = 
1.38 ± 0.02 mag for B037. Using spectro s copic m etallicities to 
predict the intrinsic colours, IBarmbv et al. I d2002bh rederived the 
redden ing value for this GC, E(B - V) = 1.30 ± 0.04 mag. Re- 
cently, IMa et al. I d2006al) determined the reddening and age of the 
B037 by comparing multicolour photometry with theoretical stel- 
lar populatio n synthesi s model s. The reddening towards B037 de- 
termined by IMa et al. I d2006al) is E(B - V) = 1.360 ± 0.013 
mag. The reddening value for B037 determined in this paper is 
E(B — V) = 1.21 ± 0.03 mag. It is clear that the consistent 
reddening values for B037 from different references co nfirm that 
this clu ster suffers from very large extinction. In fact, IMa et al.l 
(2006c) showed the dust lane across the face of the cluster using 
an HSTI ACS image, whi ch may partial ly account for its very large 
reddening value (see also lCohen 120061) . 
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